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ABSTRACT 

We study the mass-to- flux ratio M/^ of clumps and cores in simulations of 
supersonic, magnetoliydrodynamical turbulence for different initial magnetic field 
strengths. We investigate whether the M/$-ratio of core and envelope, TZ = 
{AI /^) core /(M/^) envelope Can be uscd to distinguish between theories of ambipolar 
diffusion and turbulence-regulated star formation. We analyse TZ for different Lines- 
of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average 
and median values of |7^| for different times and initial magnetic field strengths are 

typically > 1, 2) the average and median values of \TZ\ saturate at \TZ\ « 1 for smaller 
magnetic fields, 3) values of \TZ\ < 1 for small magnetic fields in the envelope are 
caused by field reversals when turbulence twists the field lines such that field com- 
ponents in different directions average out. Finally, we propose two mechanisms for 
generating values |7?.| < 1 for the weak and strong magnetic field limit in the context 
of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo 
leads to a significantly increased flux in the core and we find \Tl\ < 1. Second, in the 
strong field limit, field reversals in the envelope also lead to values \TZ\ < 1. These 
reversals are less likely to occur in the core region where the velocity field is more 
coherent and the internal velocity dispersion is typically subsonic. 



INTRODUCTION 



Understanding star formation is a fund amental problem for 
theoretical astrophysics fsee re views bv lMac Low fc KlessenI 
12004 iMcKee fc OstrikeJ 120071 ). For several years, the idea 



was that star formation is mainly reg ulated by the mag- 
netic field and ambipolar diffusion fe.g. lMouschovia!3ll984l : 
IShu et al]| 19871 ). In this model, the neutral particles, not di- 
rectly affected by the magnetic field in a gas cloud, slowly 
move inwards to the gravitational centre, while a signif- 
icant part of the magnetic flux remains in the envelope 
of the collapsing cloud. Alternatively, star formation could 
be re g ulated by supersonic t urbul ence (iMac Low fc KlessenI 
l2004l ). |Padoan fc Nordluiidl (| 19991 ) and lPadoan et all (|2004l ) 



developed a super- Alfvenic, turbulent model of dark clouds. 
Those turbulence-regulated models of star formation predict 
that clumps and cores form through turbulent compression 
at the intersection o f shocks (|Ballesteros-Paredes et al]|2003l : 
iKlessen et "ai]|2005l ). 

Both models, ambipolar diffusion and turbulence 
try to explain the relat ively low star formation rate 
(IZuckerman fc Evanj [l973 ) observed in the Galaxy, how- 
ever, the physical processes put forward in the two models 
are fundamentally different. The main problem is that the 



amount of turbulent, kinetic energy and the amount of mag- 
netic energy is typically observe d to be of the sa me order of 
magnitude in interstellar clouds (|Crutcheij|l999l ). Thus, it is 
still an open question to which extend ambipolar diffusion 
and supersonic turbulence regulate star formation. 

In order to test the two st ar formation models described 
above, ICrutcher et al.l (|2009l ) introduced a new quantity, 
TZ — {M core /{M/^) envelope, which is the mass-to-flux ra- 
tio of the core and envelope of a dense clump. If ambipolar 
diffusion plays the central role in the process of star for- 
mation, we would expect a value of \TZ\ > 1. This is be- 
cause the clouds are initially supported by the magnetic 
field, but then ambipolar diffusion finally leads to an in- 
crease of the mass relative to the flux in the centre (core) 
of the cloud. Thus, a significant fraction of the magnetic 
flux remains in the envelope during the contraction, and 
consequently \TZ\ > 1. In contrast, if the clumps form in 
a super- Alfvenic, turbulent medium, we would expect val- 
ues of \TZ\ < 1, i.e., a mass-to-flux ra tio on average beini 
higher in the envelop e than in the core (|Lunttila et al.ll200i 
ICrutcher et al.l |2009| ) . This is attributed to field reversals 
in the envelope of the cloud, due to the larger amounts 
of turbulence there, compared to the denser core. There 
are three reasons to expect \TZ\ < 1. First, the envelopes 
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of the cores typically have a muc h higher turbulent ve- 
locity dispersion t h an the inter i or jBenson fc Mverd 19891: 
Andre et all l2007l: iLada et all l2008l: iBeuther fc Hennind 



20091 : ISmith et al.ll2009l : iFederrath et al.ll201(]| '). Second, the 
magnetic field in the core is stronger than in the enve- 
lope due to the compression of the field lines. Both as- 
pects make field reversals more likely in the envelope of 
the cloud, such that magnetic field lines can cancel out 
there, leading to \TZ\ < 1. The third important mecha- 
nism leading to |7^| < 1 is small-sca le dynamo action (e.g., 
[Brandenburg fc SubramanianI [ioOSi ) in the core. This pro- 
cess leads to an efficient amplification of the magnetic flux 
in the cor e, relative to the env elope as shown in the si mu- 
lations bv lSur et al] (|2010l ') and'Fede rrath et all (l2011bl '). or 
investigated analytically bvlSchlcicl ier et al.l ( 201of K 

The idea of this paper is to perform supersonic, mag- 
netohydrodynamical (MHD) simulations of turbulence in 
dense clouds, and to analyse the relative mass-to-flux ratios 
in the core and envelope, TZ, of clumps identified in the sim- 
ulations. We use a large statistical sample in PPP (Position- 
Position-Position) and PP (Position-Position) space to ob- 
tain statistically significant results. We take into account dif- 
ferent initial magnetic field strengths (i.e., weak and strong 
magnetic fields) and investigate the dependence of our re- 
sults on the resolution of the simulations. In addition, we 
compare two different methods for computing TZ, in order 
to test the validity of our concl usions. Finally, we com pare 
our results to those given by ICrutcher et al.l (|2009h and 
iLunttila et al.l (|2008l ). and discuss the implications of our 
results on 7?. as a measure to distinguish between ambipolar 
diffusion and turbulence-regulated star formation. 

In section [2] we describe the MHD simulations, the 
clump finding and analysis, and introduce two different 
methods for computing TZ. In section [31 we show that our 
turbulence simulations produce clumps with mean values of 
|7?.| > 1 for runs with relatively strong initial magnetic fields. 
Runs with weak magnetic fields produce a larger fraction of 
cores with \TZ\ < 1, but the mean value is still slightly higher 
than unity for both analysis methods of \TZ\. In section |3] 
and[S] we present our conclusions and describe two possible 
mechanisms, field reversals and the small-scale dynamo, for 
generating values \TZ\ < 1. 



2 METHODS 

In the following, we describe our numerical methods used to 
model supersonic, magnetohydrodynamical turbulence and 
the analysis performed to define cores and magnetic field 
strengths along the Line-of-Sight (LoS), and to compute 
mass-to-flux ratios in the cores and envelopes found in the 
simulations. 



2.1 MHD simulations of driven, supersonic 
turbulence 

We computed numerical solutions of the compressible, three- 
dimensional, ideal magnet ohydrodynamical e quations with 
the grid code FLASH v2.5 (|Frvxell et al.ll2000l '). here written 



in a form where the permeability constant /io = 1: 
dtp + V ■ (pu) = 0, 

dt{pu) + V • (pu(g)u - B®B) -I- Vp* = pF, 
dtE + V ■[{E + p^)u - (B-u) B] = 0, 
dtB + V ■ (u(g)B - B(g)u) = 0, 
V ■ B = 0, 



(1) 



where p, u, = p + (1/2) |B| , B, and E = peint + 
(l/2)p|u|'^ -I- (1/2) IBI'^ denote density, velocity, pressure 
(thermal plus magnetic), magnetic field, and total energy 
density (internal, kinetic, and magnetic), respectively. The 
MHD equations were closed with a polytropic equation of 
state, p = cfp, such that the gas remains isothermal with 
constant sound speed Cs = 0.2 km s^^, assuming a con- 
stant gas temperature of 11.2 K. To drive turbulence, we 
apply the forcing term F as a source term in the momen- 
tum equation above. The forcing term is modelled with 
a stochastic Ornstein-Uhlenbeck process ( Eswaran fc Pop j 
1 19881 : ISchmidt et"al]|2009l : IFederrath et al lbOOgj ). such that 
F varies smoothly in space and time with an autocorrela- 
tion time equal to the eddy-turnover time, T = L/{2A4ca) 
at the largest scales, L/2, in the periodic simulation do- 
main of size L = 4pc. M = Urms/cs denotes the root- 
mean-squared (rms) Mach number, the ratio of the rms 
velocity and the sound speed. All models were driven to 
an rms Ma ch number of ~ 10, typical for interstellar 
clouds fe.g.. ICsengeri et al. 20 111). Turbulence is fully devel- 
oped at t = 2r (ISchmidt et al.ll2009l : IFederrath et al.ll20"lol : 
IPrice fc FederrathlbOlof T We thus analyse our results in this 
statistically steady regime for t = 2.0, 2.4, and 2.8 T to ex- 
plore the temporal variations of our results. 

The turbulent forcing is constructed in Fourier space 
such that kinetic energy is only injected at the smallest 
wave numbers, 1 < |k| L/2tt < 3, i.e., the largest scales. 
Construction in Fourier space allows us to decompose the 
force field into its solenoidal (rotational) and compressible 
(dilatational) parts. In this study, however, we only use 
solenoidal (divergence-free) forcing of the turbulence and 
leave the study of the mass-to-fiux ratios in compressively 
driven MHD turbulence for future work. 

We use d the new HLL 3R scheme for ideal MHD, de- 
veloped by IWaaganI (|2009|). and tested extensively in the 
FLASH code bv IWaagan et al.1 (|201ll ). The scheme makes 
use of a novel approximate Riemann solver for ideal MHD 
(iBouchut et al. I I2OO7I . I2OIOI ) that preserves positive states in 
highly supersonic MHD turbulence. To explore the influence 
of different initial magnetic field strengths, we performed 
simulations with an initial plasma /?o = 2po/-Bo = 0.01, 0.1, 
1, 10, and 100, corresponding to initial field strengths of 
Bo = 44, 14, 4.4, 1.4, and 0.44 /xG. AU models were evolved 
on a fixed grid with 256"^ grid zones. For the model with 
Po = 1 we additionally performed a resolution study with 
128^, 256^, and 512^ grid cells in Appendix|X| 

Figure [1] shows the time evolution of the plasma beta 
and the Alfvenic Mach number. The rms sonic Mach num- 
ber is not shown, but settles around « 10 within the first 
two turbulent turnover times, t >2T. After that, the tur- 
bulence is fully developed, and thus we restrict our analysis 
to t = 2.0, 2.4, and 2.8 T to explore the temporal variations. 
Figure [1] shows that all models are super- Alfvenic, except 
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Figure 1. Top; Evolution of the ratio of thermal to magnetic 
pressure, plasma /3, for all MHD models with different initial val- 
ues, /3o- Bottom: Same as top panel, but for the Alfvenic Mach 
number. 



for the run with fio — 0.01, 
Mach number of Ma ~ 0.8. 



which approaches an Alfvenic 



2.2 Core analysis of the numerical MHD 
simulations 

In the following we wi U use the definition of a clump and a 
core as it was given bv lBergin fc Tafallal (|2007l V They define 
a clump as an object with a size of 0.3-3 pc and a core as an 
object with a size of 0.03-0.2 pc. This means, that in every 
simulation, we will extract clumps that contain cores which 
will be used for further computations. 

We consider the full 3-dimensional information in the 
Position-Position-Position (PPP) as well as the correspond- 
ing surface density maps, i.e., the 2-dimensional projection 
of Position-Position (PP), and we note that only the latter is 
accessible to the observations. For the PPP case in our sim- 
ulation box with 256"^ grid cells, we look for the densest cell 
and select a small cube with an edge length of 15 grid cells 
around it (the density maximum is centred in the middle of 
this cube). The only information that we need for our calcu- 
lations are the density and the magnetic field components. 
For our resolution studies with 128'^ and 512^ grid cells, the 
small cubes are also scaled with a factor of 2 in each spa- 
tial direction, so that we have 7^ and 29^ grid cells for each 
clump/core, respectively. The small cube is extracted out of 
the simulation box such that it leaves an empty region in our 
box that will not be used any further. Then we find the next 



dense region and extract another clump and so on. In total, 
we select as many clumps out of each simulation box as we 
can find at any given time and initial magnetic field strength 
by their density maximum for the PPP measurements, un- 
der the assumption that the peak density does not fall below 
p — 20 p, where p is the mean density in the computational 
domain. We are able to extract about 100 clumps out of each 
simulation snapshot for the PPP case. For the PP measure- 
ments we choose a column density threshold of E = 2E, 
where E is the mean surface density in the simulation. With 
this threshold, we are able to extract about 40 clumps out of 
each simulation. The mean (column-) density of the clumps 
always ranges between 10^^^ and lO^^'' gcm^'^ for the PPP 
case and between 0.1 and 1.2 gcm^^ for the PP case. 

Our clump selection algorithm is somewhat arbi- 
trary. The selection of overdensities in a complex, fila- 
mentary structure, typical of molecular clouds, is however 
a general problem (|Smith et al.l l2008l : iPineda et al] |2009| : 
ISchmidt et~alll2010l ). In order to test at least one other way 
of defining our cl umps, we use clumpfin d , a 'friend-of-friend' 
algor ithm (e.g., IWilliams et al.l 1 19951 : iKlessen fc BurkertI 
|2000| ). and perform the same analysis. Clumpfind extracts 
clumps that are identified as connected regions. We locate 
the density maximum in each clump and define a grid cell 
being part of a core, if the density is larger than two thirds 
of the maximum density in the clump. In analogy, the en- 
velope is defined as all grid cells with density below one 
third of the maximum density. In this way, we test differ- 
ent core-to-envelope volume ratios. Using clumpfind affects 
the particular selection of clumps and hence their individual 
properties. However, we find that the overall statistical prop- 
erties of the clump ensemble are not significantly affected by 
using the clumpfind selection algorithm, and hence our main 
conclusions re main intact. 

Following ICrutcher et~aLl (|2009l '). we define the quan- 
tity TZ, which compares the mass-to-fiux ratio of core and 
envelope of our clump. 



7^ : 



[-^-^/ ^jenueZope [^/ Losi envelope 



(2) 



where the column density E denotes the Line-of-Sight (LoS) 
integral over the density, E = Jj^^g pdz, which is computed 
along any LoS of our homogeneous grid and BloS is the 
magnetic field component in the LoS-direction in which we 
observe our clumps. For a more realistic treatment of this 
scenario, all magnetic field components in each LoS are 
mass-weighted with the density of each grid cell, to make 
our computations more comparable to observed Zeeman- 
splitting measurements: 



BloS ~ 



pBzdz 



LoS 



(3) 



To calculate TZ, we first have to define a region for the core 
and the envelope, respect ively. Our aim i n gene ral is to ap- 
ply the method given by ICrutcher et all (|2009l ). who used 
one telescope beam for the core and four somewhat larger 
telescope beams for the envelope. 

Since we have many more beams (grid cells) available in 
the simulations, we use them to increase the statistical sig- 
nificance of the measurement. Figure[2]shows how we defined 
core and envelope in our clumps. After contraction (integra- 
tion) of each 3D-cube along the Line-of-Sight, we choose a 



© 0000 RAS, MNRAS 000, 000-000 



Bertram et al. 




compute the logarithmic average for A'" absolute field com- 
ponents in the envelope: 

10^E.i°g|^.l 



(4) 



Figure 2. Definition of core and envelope in our clumps for a 
resolution of 256'^ grid cells. The density peak is located in the 
centre, the blue cells define our core size, while the green cells 
define our envelope. 





PPP 


PP 


Cells in diameter 
Physical Diameter 
Mass 
Density 


5 

0.08 pc 
0.3-3 Mq 

10-i''-10-20gcm-3 


5 

0.08 pc 
2-29 Mq 
0.1-1.2gcm-2 



Table 1. Mean properties of our selected clumps for a resolution 
of 256^ grid cells for the PPP and PP case. 



circle core size of 5 grid cells in diameter around the den- 
sity maximum for the 256^ simulation (3 cells for the 128^ 
simulation and 11 cells for the 512'^ simulation) and neglect 
the complex geometry of the cloud. The envelope, which is 
treated as a thin shell i nstead of four telescope beams as 
done bv ICrutcher et al.l (|2009l ). has a size of 2 grid cells (1 
cell and 5 cells). These values refiect the average sizes of 
cores and envelopes in our simulations. Our results do not 
depend significantly on this choice of grid cells of core and 
envelope. Table[T]gives a short overview of some mean prop- 
erties of clumps for a resolution of 256"^ grid cells for the PPP 
and PP case. An alternative method for computing TZ will 
be described in the next section. 



2.3 Two different methods of computing TZ 

In general, 7?. is a statistical quantity and there are sev- 
eral possibilities to calculate its value (see, e.g., t he discus- 
sion about M/$ in lVazguez-Semadeni et "al]|201ll '). Here we 
adopt and compare two different methods. First, we com- 
pute the average of the magnetic field components in the 
LoS in the core and envelope, BloS = jj'^^Bi, and after- 
wards calculate TZ by using eq. ([2]). Second, we obtain TZ 
pixelwise, i.e., we select one pixel of our core and one pixel 
of the envelope and compute TZi only for those two pixels. 
We do that for all M possible combinations of pixels of core 
and envelope and take the logarithmic average of all abso- 
lute values afterwards, because we obtain a very wide range 
of 7?.- values, which are rather logarithmically than linearly 
distributed. We do the same for the magnetic field and also 



Our aim is to compare the results of both methods of com- 
puting TZ and BloS with each other. 

2.4 Computation of field reversals 

For our first analysis method, we also define the following 
quantity, X{N), which tells us how many field reversals ap- 
pear along the Line-Of-Sight (LoS): 

A''-!- - A'' 

Here A'^+ is the total number of cells with a positive sign of 
the magnetic field component in LoS in the envelope and 
the total number of cells with a negative sign of the magnetic 
field component in LoS in the envelope. A'^ — + N- is the 
total number of grid cells (the normalisation) counted in the 
envelope. By definition, X can only lie in a range between 
-1 to -f 1, corresponding to either cells only with a negative 
or positive magnetic field component. A value of X close to 
zero indicates that there are many field reversals in the enve- 
lope, i.e., nearly as much cells with a positive and a negative 
sign. In the following, we only consider the absolute values 
of \TZ\ and \Blos\- However, we emphasise that both quan- 
tities are computed with their individual signs (per beam 
and per cell) taken into account. This is important, because 
of possible cancellations of the B-field along the LoS due to 
field reversals. 



3 RESULTS 

3.1 PP column density map and large scale 
structures 

As an example, Figure[3]shows a logarithmic column density 
map of our simulations with /3o ~ 0.01 (left) and /3o = 100 
(right), i.e., for a very strong and very weak initial magnetic 
field. The left plot shows clear filamentary structures that 
are oriented preferably along the z-direction in which our 
initial field was established. For comparison, the right plot 
shows a rather isotropic structure, which is caused by weak 
magnetic field lines that can easily be tangled by turbulence. 
Also labelled are the positions of the density peaks and the 
maximum diameter of the envelopes of our extracted clumps 
for the PP case. 

3.2 Time evolution 

When comparing different simulation snapshots, we do not 
find any systematic time-dependence in the distribution of 
our clumps for PPP and PP for any /3o. To illustrate this 
point. Table [2] shows the statistical core values for different 
times (t = 2.0, 2.4, 2.8 T) for two extreme initial values of 
Po = 0.01 and /3o = 100 in the PPP case, computed with the 
first analysis method. On the basis of a la-error interval, no 
significant change in the distribution of our clumps (except 
for the usual statistical fluctuations) with respect to time 
are found. 
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Figure 3. Logarithmic column density map computed along the LoS in y-direction for an initial plasma beta of /9o = 0.01 (i.e. for a 
very strong field, left), and = 100 (i.e. for a very weak field, right) at t = 2.0 T. For /3o = 0.01, one can see the outstanding z-direction 
of the magnetic field, while for /3o = 100, the magnetic field is so weak that turbulence can easily tangle the magnetic field lines, such 
that the overall density structure is rather isotropic. Also labelled are the positions of 40 density-peaks (dot in the middle of each circle), 
which fulfil our threshold condition and the maximum diameter of the envelopes (circles). 



3.3 Scaling of the average values of BloS and TZ 

The plots in Figure |4] and [5] are computed with our first 
analysis method (see sec. I2.3|l . Figure |3] shows a clear trend 
between the value of \TZ\ and |-Blos| in each direction, both 
for PPP and PP. The absolute values of the magnetic field 
are smaller for smaller absolute values of TZ. For a small 
initial /3o ~ 0.01 (first row) one has very strong magnetic 
field lines in the simulation in z-direction that lead to a 
very compact distribution of clumps for every time. Here 
the magnetic field lines are so strong that turbulence is not 
able to stir the medium perpendicular to that direction. To 
quantify this effect, we show the number of field reversals 
X defined in eq. ([SJ in Figure [S] such that we can finally 
compare which clumps in Figure |4] have a certain number 
of field reversals or not. As expected, we do not find any 
field reversals for the clumps with very low /?o, observed in 
z-direction, which means that X = 1 (first row. Figure [5]). 
All clumps observed in the other directions are isotropically 
distributed. This is because of our initial direction of the 
magnetic field in z-direction at the beginning of our simu- 
lation. If we increase the initial value of plasma /3, we find 
that the average value of the magnetic field in LoS of our 
clumps decreases. In Figure [l] we a lso added the observa - 
tional results of the four clouds from ICrutcher etHI lj2009h . 
L1448CO, B217-2, L1544 and Bl, which fit into our general 
trend of increasing TZ with increasing BloS- However, the 
observed values are at the lower end of our distribution. We 
do not find any significant differences between the PPP and 
PP measurements. We also varied the number of cells of core 
and envelope, as described in section [?!2l but could not find 
any significant change in the distribution of the clumps. 

We might expect that the average values of the mean 
magnetic field of the clumps shown in Figure |4] should scale 
like BLoS,i/ BloS,] = \/A)j7Am because [3 oc B"^, where 
i and j denote simulations with different initial plasma 



p. Table [3] gives an overview of the average values of the 
magnetic field, TZ, and their standard deviations and medi- 
ans. If we consider the z-direction, we should always ob- 
tain a constant ratio of BloSj/ BloSj = VW « 3.2, if 
i and j correspond to /3o = 0.01 and 0.1, 0.1 and 1 and 
so on. Therefore, the ratio from i = 0.01 to j = 0.1 is 
BloS,o.oi/BloS,o.i = 46.2/15.8 ^ 3, which fits well to our 
theoretically predicted value of 3.2. For the other ratios, 
we get values of 1.7, 1.7 and 1.8 respectively. This discrep- 
ancy comes from the fact that the magnetic field is amplified 
by th e small-scale dynamo (jBrandenburg fc SubramanianI 
l2005h in cases of high initial /3o, i.e., /3 is a function of time, 
as is the Alfvenic Mach number (see Figure [!}. 



3.4 Effect of field reversals on TZ 

Let us now analyse the consequences of field reversals in our 
clumps on the behaviour of our statistical quantity TZ for 
our first analysis method. Figure [5] shows the corresponding 
amount of field reversals for each clump plotted in Figure U 
Besides the fact that the distribution of clumps in Figure [S] 
qualitatively moves to lower magnetic field strengths as we 
go to higher values of plasma P, we notice that the standard 
deviation of B and \TZ\ observed in z-direction is getting 
bigger for lower magnetic field strength, i.e. for a higher 
plasma /3. This is caused by the fact that a weaker field 
cannot resist as well against turbulence as strong magnetic 
fields, we therefore measure more field reversals as we go to 
higher plasma /3o for the z-direction. We also observe that 
for small magnetic fields in the LoS (independent from any 
direction) we can identify more field reversals, that means 
clumps with values of X ~ 0. 
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Figure 4. Distribution of clumps in different LoS-d i rectio ns i) for Position- Position- Position (PPP) and ii) for Position-Position (PP) 
measurements and observed cores bv lCrutcher et al ] l|2009h . From top to bottom: different values of l3o (/3o = 0.01, 0.1, 1, 10 and 100). 
Prom left to right; different time steps {t = 2.0,2.4 and 2.8 T). The initial magnetic field strength for /3o is marked with a vertical line. 
Plotted is the absolute value of TZ against the absolute value of the average of the magnetic field components for a given LoS. In general 
we observe a small value of |7^| for small magnetic field strengths, that might be caused by field reversals. The stronger the magnetic 
field lines, the higher the value of |7?.|. For PPP and PP configurations, as well as for the three different times, we get statistically the 
same distribution. 
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3.5 Comparison of two different methods of 
computing TZ 

In Section [2.31 we have described two diflterent methods of 
computing TZ for getting a statistical distribution of clumps 



in the B — 72.-scatter plot. Depending on how we average our 
values we find some differences in the scatter plots. Figure [6] 
shows the results when using our second analysis method 
for the same clumps used above. The observational mea- 
surements plotted in Figure [6] are computed according to 
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Table 2. Mean, median and standard deviation for all direc- 
tions of the magnetic field component and |7?.| for 256^ cells for 
/3o = 0.01 and /Jq = 100 for the PPP case, computed with our 
first analysis method described in section [21 From top to bottom 
(separated by a line): values for different times, t = 2.0, 2.4, 2.8 T. 
All values of B are given in /jG and time in T. 
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Table 3. Mean, median and standard deviation for all directions 
for the absolute magnetic field component and |7^| for 256^ cells 
and for t = 2.0 T for the PPP case, computed with our first analy- 
sis method described in section [2] From top to bottom (separated 
by a line): values for an initial plasma j3 of 0.01, 0.1, 1, 10 and 
100. All values of B are given in ^G. 
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Table 4. Mean, median and standard deviation for all directions 
for the magnetic field component and |7^| for 256'^ cells for /3o = 
0.01, 1, 100 for the PPP and PP case [t = 2.0 T), computed with 
our second analysis method described in section [2] All values of 
B are given in /iG. 



method 2 for the raw data eiven in ICrutcher et~all (|2009l . 
Tab. 1). Table 2] gives average values for all clumps in Fig- 
ure H For both, PPP and PP, we again find mainly |7?.| > 1 
for different /3o- As seen in Table [3] for our first analysis 
method, we also find a saturation of \TZ\ ~ 1 for high /3o 
for our second method. For example, the average PPP val- 
ues for the x-direction at different /3o are \TZ\ ~ 4.0,3.4, 1.3 
(first method) and \TL\ « 5.5, 3.6, 1.5 (second method) for 
/3o = 0.01, 1, and 100. The same trend can be observed for 
the other LoS-directions. 



4 DISCUSSION 

Here we propose two physical mechanisms to generate val- 
ues of 1 7?. I < 1 and compare to existing studies. We dis- 
tinguish between the strong and weak magnetic field limit. 
We note tha t with strong fie lds we do not refer to the clas- 
sical model (|Shu et al.lll987l ) where the energy density ex- 
ceeds all other forms of energy and ambipolar diffusion is 
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Figure 6. Same as FigureU but for the second analysis method 
describes in section [231 Only t = 2T and /3o = 0.01, 1 and 100 
are shown. 



needed to form stars. Instead, by the strong magnetic field 
limit, we refer to field strengths that are close to satura- 
tion as expected by fundamental energy equ ipartition argu- 
ments (jBrandenburg fc SubramanianI |2005| ). Detailed sim- 
ulations of the dynamo process under physical conditions, 
including magne t ic field diffusion and viscous dissipation by 
iFederrath etahl |2oili) predict the energy density of the 
field to be of order of 10 percent of the turbulent kinetic 
energy density for the typical transonic to supersonic Mach 
numbers in the interstellar medium. 



4.1 Strong magnetic field limit 

Molecular cloud cores are thought to for m at the stagna- 
tion p oints of convergent fi ows (see, e.g., iHennebelle et all 
l2008l : iBaneriee et"airi2009l ) in the turbulent interstellar 
medium, where the supersonic t urbulent moti ons domi- 
nate the cloud on large scales (|Larsonl Il98ll ) and be- 
come more coherent motions with subsonic velocity dis- 



1989: Goodman et al.l 19981: Barranco & Goodman 1998: 


Andre et al.l 2007 


: Lada et al. 2008: Beuther & Henning 


200C 


: ISmith et al.l 


20091: IFederrath et al.ll2010l: iPineda et al. 


201C 


) which is comparable to the diameter of our simulated 



clumps of 0.08 pc (see Table [T|. Because of the more coher- 
ent flow pattern, the magnetic fleld lines are less strongly 
twisted in the central core compared to the outer envelope. 
In addition, the fleld strength in the core is larger than in 
the envelope due to the compression of the field lines as the 
density increases. Both aspects make field reversals more 
likely to occur in the envelope, where magnetic field lines can 
cancel out, than in the inner core. As a result we measure 
\TZ\ < 1. In principle \TZ\ < 1 could also result from geometri- 
cal effects and observation al bias (|Vazauez-Semadeni et all 
l2005l : ICrutcher et "al]l2009l ) . 

4.2 Weal4 magnetic field limit 

Our statistical analysis of turbulent core formation in sec- 
tion [3] yielded average values of [R.\ > 1. Only for the small- 
est considered initial fleld strengths, we found a significant 
number of cores with \TZ\ < 1, but still with an average \TZ\ 
very close to unity. However, as discussed above, these calcu- 
lations focus on a regime where the magnetic field is close to 
the saturation level. It is also important to consider systems 
where the fleld strength is orders of magnitudes below this 
value. To do so, we obtained simulat ion data from lSur et al.l 
(|2010t ) and IFederrath et all (12011^ 1 where the initial fleld 
strength was only B = 10"^ G and then was amplifled by a 
factor of 10* by the small- scale turbulent dynamo drive n by 
gravitational collapse (see iKlessen fc Hennebeli3l2010l , , for 
a general discussion of the energetics of turbulence generated 
by gravitational contraction). 

Imagine a very weak initial magnetic field that is 
amplifled due to the stretching, twisting, and folding of 
magnetic fleld lines, the so-called small-scale dynamo (see 
iBrandenburg fc SubramanianI [20051 , for a review). The tur- 
bulent dynamo is a process by which turbulent kinetic en- 
ergy is converted into magnetic energy, by packing fleld 
lines closer together. This process always works, if turbu- 
lent kinetic energy is injected into the system, and as long 
as the magnetic energy is still smalle r than the kinetic en- 
ergy, from which the d ynamo feeds. ISur et al] (|2010l ) and 
IFederrath erahl (|2011bl ) showed recently that turbulent dy- 
namo amplification also works in a collapsing, magnetised 
core. Thus, even if \TZ\ > 1 initially, we expect that \TZ\ be- 
comes smaller than one, as the magnetic field gets amplified 
in the core. Figure [7] shows the time evo lution of \ TZ\ for 
the turbulent, magnetised c ore studied in ISur et all (|2010D 
and IFederrath et al.l (|2011bh . The core starts off with a very 
weak initial magnetic field and \TZ\ — 5.8 at time t = 0. The 
initial turbulence prevents the core from collapsing imme- 
diately, and makes \TZ\ decrease to around unity at r « 4, 
when core collapse sets in. This initial decrease of \TZ\ is 
due to field reversals as explained for the turbulent mod- 
els in section 13.41 During the collapse phase, however, the 
magnetic field is amplified further, and for r > 4, we find 
values \TZ\ < 1, as expected, due to the dynamo-amplified 
flux in the centre of the core. We conclude that magnetic 
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Figure 7. Time evolution of \'R,\ for the t urbulent, weakly mag- 
netised, collapsing Bonnor-Ebert sphere of lSur et al. I bOlOl) . The 
initial turbulent phase, r < 4, leads to a decrease of |7^| to 
around unity, duo to field reversals. The magnetic field gets fur- 
ther a mplified by the gravity- driven, small-scale turbulent dy- 
namo jFederrath et al.l [ioiThl) during the subsequent collapse 
phase, T > 4, which leads to |7?.| < 1. 



fie ld amplification durin g si'avitational collapse, as reported 
in lSur et all (|201(]| ') and iFederrath et all (|20inj ). can be an 
important mechanism to explain the low values o f |7^| < 1 
found in the core observations bv lCrutcher et all (|2009l '). 



4.3 Relation to existing studies in the literature 

On avera ge, our results s eem to show the opposite trend 
found bv iLunttila et al] (|2008t ). who observed that most 
clumps have values of \TZ\ < 1 for larger magnetic fields 
and values of \TZ\ > 1 for lower magnetic fields. However, we 
find approximately 70-90% of our clumps for /So = 0.01 and 
40-60% of our clumps for /?o = 100 having values of |7?.] > 1, 
despite of a very broad Icr-interval. Thus, for small magnetic 
field strengths, we find more clumps with small \TZ\ than for 
higher field strengths. We interpret this as being a result of 
field reversals, which naturally leads to smaller values of |7?.|, 
if the field is more easily tangled. In addition to that, small 
fields can be amplified by the small-scale turbulent dynamo 
action, again leading to a predominately smaller values of 
\TZ\. Both aspects lead to the trend we observe in our statis- 
tical analysis, i.e., that \TZ\ is smaller for smaller magnetic 
fields. 

This trend is basically in agreement with the 
ICrutcher et~all (|2009l ) cores although only four clumps are 
analysed. Any trend seen in their observations might thus be 
due to low-number statistics and does not necessarily reflect 
average clump properties. Similar holds for the average val- 
ues. The average v alues of |7^| in our sim ulations are mostly 
larger than 1, but ICrutcher et al] (|2009l ) only observes val- 
ues \TZ\ < 1, with some dependency on the analysis method 
(see sec. 13. 5p . The question remains whether this is caused 
by a physical or statistical effect, because of lack of observa- 
tional data. In fact, our results show that the quantity \TZ\ 
is not necessarily a good statistical measure to distinguish 
between ambipolar diffusion theory and turbulent theory of 
star formation, because of the large distribution of \TZ\ with 
values of \TZ\ both larger and smaller than unity. In fact. 



we find that most of our clumps have values with \TZ\ > 1, 
even in a turbulent environment without ambipolar diffusion 
acting in our simulations. 

One might argue that numerical diffusion could have 
had a similar effect as ambipolar diffusion, and that the 
fact that we obtain mean values |7?.| > 1 is a numerical ef- 
fect. This, however, can be safely excluded as shown in Ap- 
pendix |X] There, we show that the mean properties of our 
clumps do not show a systematic de pendence on resolution. 
In addition. IFederrath et al] (|2011bl . App. C) explicitly mea- 
sured the numerical diffusion of the present MHD scheme. 
We thus find that numerical effects have no impact on our 
general conclusions. 

A difficulty of interpreting these results is the pro- 
jection of clumps along the LoS-axis for the PP case. 
The projection leads to blending of structures in the LoS, 
so that a clump found in the column density map (PP) 
could also be a superposition of more than one object (or 
filament) that are spatially separate d from one another 
(IBallesteros- Paredes fc Mac LowllioO^ '). This is the reason 
why we also analysed the PPP case and compared both. 

Finally, we would like to mention the observa- 
tional challenges of measuring BloS in the envelope 
of a clu mp and the recent discussions in this context 



Mouschovias fc Tassid [20091 : ICrutcher et~al]|2010l : 



iMouschovias fc Tassisll2010l l 



5 SUMMARY AND CONCLUSIONS 

We have analysed MHD simulations of driven, supersonic 
turbulence with different initial magnetic field strengths 
(^0 = 0.01, 0.1, 1, 10, and 100) at different tur- 
bulent turnover times, t — 2.0, 2.4, 2.8 T. Our aim 
was to analyse the statistics of the parameter TZ — 
{M /^) core/ {M /^) envelope for a Statistically significant num- 
ber of clumps. To compare our predicted values with obser- 
vational data, we distinguished between the PPP and PP 
case. For each model, we extracted cores that were found by 
their density and column density peaks and calculated \TZ\ 
and \Blos\ for each of them. We also introduced a quantity 
X{N), as a measure of the number of field reversals in the 
envelope of a clump and compared two different methods for 
computing TZ. 

We identify two primary physical processes by which 
TZ <1 can be achieved in turbulent magnetised clouds: 

• In the strong magnetic field limit, field reversals are 
more likely in the envelope than in the core (see sec. 14. ip . 
such that magnetic field lines can cancel out there, leading 
to \TZ\ < 1. 

• In the weak magnetic field limit, the small-scale dynamo 
amplification results in an increase of the initial magnetic 
field strength in the core, thus also leading to values \TZ\ < 1 
(see sec. l4.2p . 

In addition, we report the following findings: 

• We saw no significant time evolution in our distribu- 
tions on the basis of la-interval for all /3o and LoS-directions. 
This is because in all our simulations, turbulence was already 
fully developed and the magnetic field was saturated. 

• We did not measure any significant differences between 
PPP and PP, neither for any LoS nor for any /3o- We found 
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that those values with \TZ\ < 1 are mainly generated by 
field reversals in the clumps, which lead to a higher average 
mass-to-flux ratio in the envelope than in the core. 

• Our distribution has average values of \ 7V\ > 1, but has 
a larg e standard deviation in \R.\ and |i?r,n,c;|. ICrutcher et alj 
12009) observes average values of \TZ\ < 1, for four different 
cores/clumps. However, their four cores are consistent with 
our basic trend, but are located at the lower end of our 
distribution. 

• We find a simila r trend as seen in the observations by 
ICrutcher et al.l (|2009l ). namely that \TZ\ is smaller for smaller 
magnetic field strengths. As for the mean values of the ob- 
served [R\, however, the trend does not necessarily refiect 
the average properties of typical clumps. 

• The actual values of \TZ\ depend slightly on the analysis 
method. Our second method (sec. 12. 3p gives slightly larger 
values of \TZ\ on average. Applying the s e cond method to the 
four clumps observed in lCrutcher et all (|2009l ). we find that 
one of the four clumps has \TZ\ > 1, compared to the first 
method, for which all four clumps have \TZ\ < 1. 

For further investigations and to get a better under- 
standing of the distribution of \TZ\ in terms of the aver- 
age magnetic field strength, we need additional efforts in 
both observations and simulations. While one has to observe 
many more clumps to get a better, statistically significant 
distribution, future numerical simulations must eventually 
incorporate turbulence and ambipolar diffusion, as well as 
an accurate treatment of the cores' chemical evolution to- 
gether with a proper model for Zeeman splitting and radia- 
tive transfer. 
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APPENDIX A: RESOLUTION STUDY 

In this section we present the results from our resolution 
studies (the values are computed with our first analysis 
method, sec. 12. 3p . Therefore, as an example, we focused on 
the simulation with an initial plasma /3o = 1 at f = 2.0 T for 
resolutions of 128^, 256^, and 512^ grid cells (the physical 
properties of our numerical experiment are not changed). 
Figure lAll shows that there is no significant change in 
the distribution of our clumps, except the scattering due 
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Table Al. Mean, median and standard deviation for all direc- 
tions for the magnetic field component and |7^| for a resolution of 
128'^, 256"^ and 512"^ cells (from top to bottom separated by a line) 
for the PPP case {t = 2.0 T) computed with our first method. All 
values of B are given in fiG. 

to statistical fiuctuations. Table lAll gives an overview of 
some statistical moments (mean, median and standard de- 
viation) of the magnetic field component and \TZ\ in each 
direction. The values do not differ significantly and are 
within the Icr-error-range, for example, |-Blos| ~ 8.1 ± 5.2, 
9.5 ± 6.5, and 8.8 ± 5.8/iG for the different resolutions in 
z-direction. Also for the y-direction, where the variation is 
the strongest compared to the other LoS-directions, we have 
\Blos\ = 6.9 ± 5.2, 9.3 ± 6.6, and 7.4 ± 5.8^G. Even in this 
case, the differences from the average values are smaller com- 
pared to the standard deviations of the distributions. Also 
the change of the average values of B and \TZ\ with time (as 
seen in section [3| is at most as big as the change caused 
by effects of resolution. Hence, our results do not depend 
significantly on resolution. 
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